hepth@xxx/9408057 
LA-UR-94-2075 



EXACT, E = 0, SOLUTIONS FOR GENERAL 
POWER-LAW POTENTIALS. 
I. CLASSICAL ORBITS 

Jamil Daboulf] 

Physics Department, Ben Gurion University of the Negev 

Beer Sheva, Israel 

and 

Michael Martin Nietog 
Theoretical Division, Los Alamos National Laboratory 
University of California 
Los Alamos, New Mexico 87545, U.S.A. 

ABSTRACT 

For zero energy, E = 0, we derive exact, classical solutions for all 
power-law potentials, V(r) = —^ijr v , with 7 > and —00 < v < 00. 
When the angular momentum is non-zero, these solutions lead to the 
orbits p(t) = [cos //(¥>(*) - ^> (*))] 1/m , for all n = v/2 - 1 / 0. When 
v > 2, the orbits are bound and go through the origin. This leads to 
discrete discontinuities in the functional dependence of (p(t) and (po(t), 
as functions of t, as the orbits pass through the origin. We describe a 
procedure to connect different analytic solutions for successive orbits at 
the origin. We calculate the periods and precessions of these bound orbits, 
and graph a number of specific examples. Also, we explain why they all 
must violate the virial theorem. The unbound orbits are also discussed in 
detail. This includes the unusual orbits which have finite travel times to 
infinity and also the special v = 2 case. 
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1 Introduction 



Since the birth of quantum mechanics, studying the connections between classical 
and quantum physics has been an enormous industry jE], p], ||. Both the intuitive 
and also the analytic aspects have been studied. In particular, one can observe a 
kind of "Folk Theorem." Problems that are amenable to well-defined, exact, analytic 
solutions tend to be solvable in both the classical and the quantum regimes. 

For example, the fundamental three-dimensional problems of quantum mechanics, 
the harmonic oscillator and the hydrogen atom, are those classical problems which, 
by Bertrand's theorem 0J, have exact, closed orbits for all energies. The standard, 
one-dimensional, quantum-mechanical potential problems that are exactly solvable, 
such as the Morse ||, Rosen-Morse ||, and Poschl- Teller potentials, are also 
solvable classically || |], [L(J . Further, it has even been shown that quantum problems 
which are solvable, but have certain "mathematical diseases," manifest these diseases 
already in the classical problem [11]. 



Usually, when one solves a potential system, one takes a particular potential from 
a family and solves it for all values of the energy, E. Here, we are going to do the 
opposite. We will consider the entire class of power-law potentials, parametrized as 

V(r) = - ^ = --^ } 7 >0, -oo<z/<oc, (1) 

and exactly solve it for all v with the particular energy E — 0. Note that it will be 
useful to switch back and forth between the variables v and /i, related by 

H={y-2)/2, u = 2(fi+l). (2) 

The potentials (I) are attractive for v > and repulsive for v < 0. For v = 0, the 
potential is a constant, V(r) = —7, so the particle is force-free. 

The above system is exactly solvable in both the classical and quantum-mechanical 
cases, and there are similarities in the properties of the solutions. In the present paper, 
hereafter called I, we shall discuss the very unusual and enlightening classical solutions 



for the orbits. In the following paper [12|], called II, we will solve the quantum problem 
and discuss the wave functions. 

In Appendix A we give, for completeness, the simple solutions for zero angular 
momentum. The main body of the paper concentrates on the more interesting general 
solutions, with non-zero, angular momentum. 

We begin, in Sec. 2, by defining our "dynamical" units, which allow us to simplify 
the equations and their solutions. Then we show that for v 7^ 2, the classical solutions 
have only one finite, non-zero turning radius, which we will denote by a. This a has 
opposite interpretations, depending on the value of the power index v. For v > 2 the 
solutions yield bound orbits with < r < a, but for v < 2 they yield unbound orbits, 
with r > a. In Sec. 3 we derive the general E = solutions, using two methods. In 
Sec. 4 we investigate the general properties of the bound orbits, v > 2 or fi > 0, such 
as their precessions and periods. 

Continuing, in Sec. 5 we graphically depict examples of the bound orbits, thereby 
demonstrating the general properties derived previously. Beginning with the case 
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v — 4, which is a boundary case that does not precess, we find that the orbit describes 
a circle which goes through the origin and continually repeats itself. For v > 4, the 
typical solution looks like a flower. A single orbit starts at the origin, goes out to a 
and returns to the origin. This first orbit describes what we will call a "petal." The 
following orbits then go on to describe further petals which in general are precessed 
from each other. After a number of petals, the trajectory closes if v is a rational 
fraction. The petals become narrower and narrower as v — > oo. Contrariwise, for 
4 > v > 2, the orbits become tighter and tighter spirals in and out of the origin as 
v — > 2. Once again, for v being a rational fraction, the orbits eventually close on 
themselves. 

In Sees. 6 and 7 we discuss the classically unbound orbits. These are given 
by all v < 2 as well as the special case v — 2. We close with comments on the 
classical problem and defer comparisons between the classical and quantum problems 
to paper II. Appendix B discusses the scaling properties of the solutions. Appendix 
C investigates violations of the virial theorem. 



2 Classical Units and E = Orbits 

Power-law potentials do not have a built-in length scale. To specify a power potential 
completely, two units are needed: a unit of energy, £ , and a unit of length, a. Instead 
of the unit of energy, one can use a unit of angular momentum, L , which can be 
related to £ Q , for example, by 

S = L 2 /2ma 2 . (3) 
Thus, we can parametrize the power potentials as follows: 

Vvr) = = —to — = — — , p = r a . (4) 

v ; V 2ma 2 p" ' H ' y ' 



Therefore, the dimensional coupling constant, 7, is related to the dimensionless cou- 
pling constant, g 2 , by 

7-2^-2 



Dimensionless units, such as g 2 , are convenient in quantum mechanics. There one 
has other physical constants, such as Planck's constant, H, to make a convenient choice 
of L and £ . For example, in Sec. 2 of II we choose L = h to obtain £ = h 2 /2ma 2 
as the unit of energy, but we leave a arbitrary. (We also could have chosen a to be 
equal to h/mc = A c , the (angular) Compton wave-length, where c is the velocity of 
light. Then we would have obtained £0 = mc 2 /2.) 

However, for classical orbits, it is often convenient to use dynamical units. By 
these we mean units which depend on the orbit itself and are determined by the 
initial conditions. For our purposes it is convenient to choose the dynamical units 
such that a is the turning point and L is the corresponding angular momentum for 
the E — solution. The possibility of such a choice is based on the following result: 
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(a) Let u a" be a prescribed, non-zero length scale. Then every power- 
law potential, V = — r y/r l/ with 7 > and v 7^ 2, will have an E = 
classical orbit with r = a as a turning radius. We call this orbit the 
"standard orbit. " 

(b) This turning radius corresponds to the maximum (minimum) dis- 
tance from the center of the force for v > 2 (v < 2), respectively. 

(For v = 2 the orbit has no turning point, except for the origin.) 



Part (a) of the above result follows from energy conservation: 



T + V(r) 



m ( dr\ L 2 
~2\dt) + 2mr 2 



-JL = ^(±) +U(L,r) 
r v 2\dt v ' 1 



E 



(6) 



where L 2 /2mr 2 is the centripetal-barrier potential and 



U(L,r) 



^ 7 
2mr 2 r v 



(7) 



is the effective potential. Since a turning radius corresponds to dr/dt = 0, the turning 
radius a must be a zero of U(r) for E = 0. Thus, by combining Eqs. @ and (^), we 
obtain 



U(L,r = a) 



1 



2ma 2 a u 2ma 2 \Lq 



L 2 2 
T2 -9 



. 



(8) 



Therefore, by choosing the angular momentum to be 

L = gL , (9) 

we obtain the desired "standard orbit." 

By using the dynamical unit L, instead of Lq, the potential U is simplified to 



U(L,r) 



L 2 / 1_ V 

2ma? \p 2 p v 



r 

P = ~ , 

a 



(10) 



so that 



V(r) = - — 



L 2 1 



2„v-2 



L 2 a 



2ma 2 p v 



2m r h 



[ID 
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Finally, from Eq. ((TTJ) we can read off the relation between L and the turning 
radius a: 



a 



In particular, for L = there is no turning point, except for a = oo. This special 
case is treated in Appendix A. 

Part (b) of the result follows by substituting E = into Eq. (|6]) to yield 

m ( dr\ 2 , . 

^ = "2 * £ °- < 13) 



Combining this with Eq. (|10D means that 

p 1 " 2 < 1 • (14) 

This condition determines the physical domain of the distance r: 

r < a, z/>2, (15) 
r > a , < 2 . 

Therefore, a is the apogee length for bound orbits {y > 2) and the distance of closest 
approach for unbound orbits {v < 2). 

For i/ = 2 the above arguments break down, since the effective potential U = 
(L 2 /2m — 7)r -2 does not allow a finite, non-zero turning point when E = 0. This 
special case will be discussed further in Sec. 6.1. 

Finally, we note that, having established the existence of a "standard orbit," it 
follows by scaling arguments that there are an infinity of similar orbits whose turning 
radii and angular momenta are related to those of the standard orbit as follows: 

a x = Xa , L\ = X l - u/2 L = X^L . (16) 

These scaling relations are, of course, just special cases of the more general E ^ 
scaling law. (See Refs. O, [HJ and Appendix B.) 



3 Classical Solutions 
3.1 From the orbit equation 

A direct, mathematically-inspired solution can be obtained by integrating the orbit 
equation |TR |TB[ 

r r" 1 dr 

Lf - Lf = / . (17) 

' ' ' '%(E-V)r*-1 
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The general solution has the constant ipo in it. We presciently set <^ = when p = 1 
(the turning-point condition), so that the turning point is along the positive x-axis. 

Then, setting E = makes the integral (O) doable for all v. Changing variables 
successively to p = r / a, y = p^, and x = cos -1 y means that Eq. (fl7l) can be written 

as 

fp pfc- 1 * dp _ x fv dy /-cos- 1 j, 

Lp = l7T^ = (I 17^ = -" i dx - (18) 

Therefore, the solution is 

p^ = y = cos (— /i<p) = cos (/i<p) , (19) 

or 

2 

" 2 = [cos(^)]^. (20) 

We will discuss the allowed angular variations of <p in the separate sections on bound 
and unbound orbits. 



cos 



{u-2)<p 



3.2 From the energy conservation equation 

We now give a more intuitive derivation of the solution (]2"0"D. By substituting the 
angular-momentum conservation condition 

<p = L/{mr 2 ) (21) 

into the energy conservation condition 



m 



(22) 



one obtains [0, [H| 



5) + r! 



2m(£ - V)r 4 

T 2 



(23) 



This is essentially a first-order differential equation, which can be formally integrated 
to yield the angular equation (|TT| ) of the last subsection. 

However, for E = 0, it is much more efficient to solve Eq.(P3"|) directly. Converting 
to the dimensionless variable p = r/a and substituting V from Eq. (|IlD into Eq. 
(|23T), we obtain 



dp^ 



P 



(2-2/1) 



(24) 
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For v = 4 the right-hand side of this equation is unity, so the solution is a cosine. 
This is the circular orbit p = cos ip which we will discuss in detail in the next section. 
Guided by this and the substitution y = p^ of the last subsection, we multiply Eq. 
by p 2 »- 2 to yield 

Now p M satisfies the differential equation for the trigonometric functions. Therefore, 
the general solution of Eq. ( p5|) is given by 



p M = cos p(ip — ipo) = cos 



v-2 



(26) 



or 



p = [cos n((p- (p )] 



COS [ ^— — (<£ - (y9 ) 



(27) 



The phase, <po> is the integration constant. It depends on the initial conditions and, 
as before, we will set <p to zero for our first-orbit turning-point condition. 



4 Properties of the Bound Trajectories, 2 < v or 

1 < fJL 

In Fig. 1 we show a typical effective potential U(r) in this regime. One sees that 
the value of U(r) starts from — oo at r = 0, rises through zero at r = a, reaches a 
maximum, and then decreases to zero at r — *• oo. Therefore, the E = solution for 
r = r(t) can only vary between r = and r = a. We shall now discuss the dependence 
of r on the azimuthal angle, (p. 



4.1 The first orbit 

The radius, p, is nonnegative. This means the range of the angle variable for the first 
orbit is restricted. Since, by convention, the first orbit has <p = <Po = 0' this means 
ip 1 = ip satisfies 

= <m l < = — . 28 

2/i z/-2 _ ^ _ z/-2 2/i y 1 



The corresponding curve p = p{tp l ) in the x-y plane, begins at p = for ip 1 = 
—n/2[i = —tt/(u — 2), evolves counterclockwise to p = 1 at ip 1 = 0, and then contin- 
ues on back to p = at ip 1 — ir/ (2/i) = 7r/(v — 2). The curve evolves counterclockwise 
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because we have chosen, by convention, that the angular momentum is in the positive 
z direction. (Of course, there is also a mathematically-rotated clockwise solution.) 

Every such closed circuit, beginning and ending at the origin, we will call an orbit. 
Because of the shapes of the orbits (which we will present in the next Section), when 
v > 4 we will also call an orbit a petal (of a flower) and when v < 4 we will also call 
an orbit a spiral, which actually is a double spiral since it spirals out and then spirals 
in. These two classes of orbits meet at v — 4, which is the circle that goes through the 
origin. As we will see below, the entire physical trajectory will consist of a pattern of 
either (i) a finite integral number of orbits which then repeat over themselves, when 
fi is a rational number, or else (ii) an infinite number of nonrepeating orbits, if [i is an 
irrational number. (In Sec. 5 we will discuss the classically unbound orbits described 
by v < 2 or /i < 1.) 

4.2 Solutions for later bound orbits 
4.2.1 The angles and the phase shifts 

To obtain the entire physical solution, we must connect different orbits. That is, the 
first orbit evolves into the second orbit... evolves into the kth orbit. Because the orbits 
go through the origin, both <p and ipo are discontinuous, as we will demonstrate below. 
Therefore, we will label (p and y? for the kth orbit as (p h and </?q, respectively. 

The successive orbits must be connected at the origin, r = 0, in such a way that 
the directions of the linear momentum p and the angular momentum L are continuous 
functions of time. We defined the first orbit to have its apogee, r = a, at ip 1 = 0, 
which corresponded to the choice <pl = 0. In order to obtain a physical solution, we 
must connect the first orbit to the second orbit in such a way that the tangent vector 
to the combined trajectory is continuous at the origin. 

Before the particle starts the second orbit, the trajectory goes though the origin. 
This creates a singular transition in the polar coordinates. (Since the trajectory passes 
through the origin and because p must be continuous, the direction r of the position 
vector must change sign at r = 0.) The angle advances by ir (instead of — 7r) in going 
through the origin because, in order to conserve angular momentum, the position 
vector must continue rotating in the counterclockwise direction. Then the particle 
begins its second orbit, travelling another angular distances 7r//x. Before beginning its 
third orbit, the particle goes through the origin again, advancing another n radians. 
Therefore, the angular variation of the kth orbit, ip k , is given by 

<<Ln - + (*-!))*<*>*< (^y^ + (* - 1)) tt = V k max (29) 

and one has the condition 

</4tn = <fL* + * • (30) 

However, the phase shift, (p%, must also change with each orbit. Recall, in sum- 
mary, that the k-th orbit is described in polar coordinates by 

p fc =[cos/V-^)] 1//1 , (31) 
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with 

\<P k -<Po\ <vr/2/i (32) 



or, equivalent ly, 



<fLn = - 2jl < V k < + 2jl = ¥>Lx ■ (33) 



Then, from Eqs. fl2Tf ) and (§3|), the phase shift, y?Q, is 



ifl = ^ + = (Jfe - 1) f — ^— ) IT = (k — 1) | 1 + — J 7T . (34) 



v — 2 / V // 



Mathematically this change of phase shift, ip^, with each orbit is because the orbit 



must be symmetric about any apsidal vector [19]. This is true even though the vector 
is zero length from the origin. (See the interesting application of these concepts in 
the v = 6 orbit discussion of the next section.) Further, the other apsidal vectors, 
Yk{ip k — <fo), are the symmetry axes of the kth orbits. 

4.2.2 Angle vs. phase-shift differences 

It is instructive to observe that, instead of using the label, k, for each orbit, formally 
we can write the general solution for all times t > in the form 

pit) = [ cos n(<p(t) - <f (t)) = [ coanxV) ] 1/At , (35) 

with 

X {t) = (p(t) - <f (t) . (36) 

Here, both the azimuthal angle, <p, and also the phase shift, <p , are regarded as 
functions of time. These angles have discrete jumps at the crossing times 

t k = (k — l/2)r , (37) 

where the are the times at which the particle passes through the origin and r is 
the period of one orbit. Using our same physical convention, which here is r(0) = a 
at t = with ip(0) = v?o(0) = 0, we have that the jumps are 

<p(t h + e) = <p(t k - e) + 7T (38) 

and 

ipo (t)=n(l + -)f^Q(t-t k ) , (39) 
where 0(i) is the Heaviside step function. 
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We see that both and y?o are monotonically increasing functions of time. In 
contrast, their difference, x{t)> * s a periodic function of t. For every new orbit, x 
starts with the value — 7r/(2/x) and then increases monotonically to ir/(2/i), due to 
the continuous increase in (p(t) as the particle goes through one orbit (petal or spiral). 
Next, x decreases stepwise by n/fj, at the crossing time. Therefore, as in our previous 
representation, cos fix starts and ends each orbit at p = 0. 

4.3 Precession of the orbits 

Because of angular-momentum conservation, the azimuthal angle, (p, is a monotonic 
function of time, t. If there is no precession, (p increases exactly by 2tt after one 
period, r. However, if after one period the axis of the orbit has rotated forwards or 
backwards (i.e. clockwise or counterclockwise) from 2ir by an angle P u , then this is 
the precession per period. Specifically, 

<p(t + r) - <p(t) = 2vr + P v , (40) 

where r is the period of one orbit. If we choose t to be the time at which the particle 
was at the apogee of the kth orbit, then t + r will be the time at which the particle 
reaches the apogee of the (k + l)th orbit. To do this, the particle must first rotate by 
7r/2/x to reach the origin, have its angle <p jump by tt at the origin, and finally rotate 
by another tt/2/j, to reach the new apogee. Thus, from Eq . ((40|) we have 



*< + r)-*t) = £+»+i=g + l).r = ar + ft, (41) 



so that 




P u =(--l)n=(LjQ7r. (42) 

In Fig. 2 we plot P v as a function of v. The precession is infinite at v — 2, decreases 
through zero at v — 4, and asymptotes to — tc as v goes to infinity. Therefore, the 
precession will be zero for v = 4, negative (counter-clockwise) for v > 4, and positive 
(clockwise) for v < 4. 

This also demonstrates that, for v equal to a rational fraction, the orbits will 
eventually close on themselves, and repeat. This happens when kP v is an integer 
times 2tt. 

4.4 Classical period 

The classical period, r, can be obtained by integrating the angular-velocity equation 

d i = —^—*> («) 

at mr z Top 
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where r is a convenient unit of time, 



ma 



r 



L 



(44) 



One then has that the classical period in units of r is 



rp T f 

TO 



\v-2\ 



p 2 (tp) dip 



2 1 M 



2R 



[cos(|/x|<^)]^ dip 



(45) 



(In Eq. (|45|) we have inserted absolute values around \i in the arguments of the 
integral and integrand. This allows us to include the negative \i case that we will 
return to at the end of this discussion.) 



Changing variables to x = \fj,\ip allows the integral to be evaluated as |20 



[ 2 [cosx](m) dx = 7^5(1/2, b) 



7T 



r(6) 



|/i| r(6 + i/2) 



b > 



(46) 
(47) 



where 



1 1 

Jl + 2 



2z/-4 



> 



(4* 



B(b,c) = r(6)r(c)/r(6 + c) is the beta function and we used T(l/2) = y/n. We see 
that T v is finite and well-defined for fi > or v > 2, as would be expected for bound 
orbits. 

In Fig. 3 we plot T v as a decreasing function of v. From Eq. (f4"7j ) special cases 

are 



limTl 



2+e 



(49) 



T* = — 



3tt 
T 



T 4 



7i 

2 ' 



lim Tj, 

!/— >00 



2tt 

V 
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The first equality in Eq. ( pE9|) is obtained by using |2]J the relation T(z + 1/2) ~ 
T{z)y/z ) which holds for \z\ — > oo. 

In passing, note that Eqs. ( f4T|) and (|48| ) also tell us that T), is finite and well- 
defined for the unbound, infinite orbits corresponding to fi < —2 or v < —2. This 
means that it takes a finite time for the particle to travel in from infinity, reach the 
turning point r = a, and then go back to infinity 

4.5 Geometric approximation to the area of the orbit 

Kepler's 2nd Law states that equal areas are swept out in equal times. Although 
usually considered for the gravitational potential, it actually holds for any central 
potential. The proportionality constant follows from the angular-momentum conser- 
vation equation (pl|): 



dA = -r 2 dp = -r 2 (p dt= — dt , (50) 
2 2 2m 



which is related to the angular-velocity equation fl4"3|) in an obvious way. Integrating 
this equation over a complete period, we get 



A, = ±r„ , (51) 



where once again the subscript v labels the potential we are considering. We now 
introduce a dimensionless area, A v , as the area of an orbit in units of a 2 : 



K = J ^~ 2^ T " ~2V =l T - (52) 



Thus, the reduced area, A u , is equal to one-half the reduced period, T v . 

For these v > 2 or fi > cases, we can give amusing, geometric, upper and lower 
bounds to the areas of the orbits. Recalling that the opening angle of a petal or spiral 
is 

». = - = -^, (53) 
H v — 2 

an upper bound is defined by the area of section of a unit circle with opening angle 
<& v . This bounds the real 



A v < \<$> v . (54) 
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Of course, this sectional area is equal to half the opening angle itself, measured in 
radians. 

We can also obtain another bound, by considering the quantity 

p((p) = cos pip (55) 
in comparison to the real orbital radius 

p = [cospp} 1/fl . (56) 
By using p(ip) in place of p in the integral for A u , 



A » = o /„ 2 " P 2 (<P) d ¥= n \T [cospp]^ dip , (57) 
one obtains the quantity 



= \ / 2 ; p 2 (ip) dip = $,/4 . (5? 



yly is an upper bound of the reduced area when 2<z/<4or0</i<l, because then 
p(<p) will be outside the area of the orbit. However, for 4 < v or 1 < p, A v becomes 
a lower bound, because then p((p) is inside the orbit. 

From the above, we can make special case comparisons of these bounds, similar to 
the results of Eq. ( pE9|) . Differing by a factor of 2 because of the reduced area definition, 
and written in the form "lower bound" < "exact reduced area" < "upper bound," 
one has 



\imA 2+e 

e— >0 



7TV/2 7T 1 

2i) < S = 4* 2+ «' 



(59) 



37T 71 1 

J 8 2 4 



$4 



71 



71 



A A 



< A 6 



7T L 71 1 

- = -$ 4 < - = -$ 4 

4 4 2 2 



1 7T 1 

- < - = -$f 

2 4 2 



7T 



< lim A, 



2(i/ — 2) "->«> 



7T 7T 1 
< = 

i/ i/-2 2 



The upper bound ^ v is very good as v — > oo, because the area of the orbit is 
a very thin petal. As z/ becomes smaller, the approximation becomes worse. By the 
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time v = 4 or /x = 1 this upper bound is off by a factor of 2. This is because the 
reduced area is that of a circle of radius 1/2, whereas this upper bound is equal to 
half the area of a unit circle. Both upper bounds become useless as v — > 2. This is 
because the upper bounds $t,/4 and $ v /2 are describing areas with angular widths 
approaching infinity and radius 1, whereas the true area has a radius that is spiraling 
in to zero as the angular position becomes large. 

Finally, note that the lower bound goes from being off by a factor of 2 as v — > oo 
to equality at v — 4. 



5 Examples of Classical Bound Trajectories 

5.1 Petals: v > 4 or fi > 1 

From the previous section we see that, as v approaches infinity, the orbit is an ever- 
narrowing petal of width 

= - = ^ • 60 

\i v — 2 

Similarly, the counterclockwise precession reaches — it per orbit as v —>■ oo. As v 
becomes smaller, the petals increase in width. We illustrate this with several illumi- 
nating examples. 

v = 8, three petals: We begin with the case v = 8. Here a petal is tt/3 wide and 
the precession per orbit is — 2n/3. Thus, there are three orbits before the trajectory 
closes. Note that the three petals in a closed trajectory cover only half of the opening 
angle from the origin. We show this in Figure 4. 

v = 7, ten petals: For v = 7 a petal is 2tt/5 wide and the precession per orbit is 
— 37r/5. This means that, before the trajectory closes, there must be ten orbits and 
the precession goes around three times. 

v = 6, perpendicular lemniscates: The case v = 6 is very interesting. The width 
of a petal is n/2 and the precession is — 7r/2 per orbit. Here, the width of a petal and 
the precession are exactly such that there is no overlap and also no "empty angles." 
It takes four orbits to close a trajectory. This is shown in Figure 5. We see that 
the physical solution consists of two perpendicular lemniscates (figure-eight curves 
composed of two opposite petals). This conclusion is in contradiction to one given 
in the literature [^| . There, a single lemniscate was predicted. From the physical 
arguments we have given above, this is not valid. This incorrect solution resulted 
from studying the equation for the square of the orbit. 

v = 5, six overlapping petals: The petals are 2tt/3 wide and they precess by — n/3 
per orbit. This means there are six petals in a closed trajectory (flower). Note that 
here a petal is still so wide and the precession is so small, that the successive petals 
overlap. 
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5.2 Circle through the origin: v = 4 or ji = 1 



For z/ = 4, the solution is well known [23|. It is a circle that starts at the origin, travels 
symmetrically about the positive x-axis, and returns to the origin. The precession is 
zero, so the orbit continually repeats itself. In Fig. 6 we show the orbit. 

5.3 Spirals: 2<z/<4or0<^<l 

As v becomes less than 4, we can think of a petal obtaining a width greater than ir, 
i.e., an orbit consists of two spirals, one out and one in, at opposite ends of the orbit. 
As v approaches 2. the spirals become tighter and tighter and the precession (now 
clockwise) becomes larger. In fact, the spirals' angular variation as well as the orbit's 
precession both become infinite in magnitude as v approaches 2. 

Consider the special case v = 3. The width of the double-spiral orbit is still given 
by the formula for <& u , and is 2n. Therefore, the first orbit begins and ends towards 
the negative x-axis. The precession is ir, so the trajectory closes after two orbits. We 
show this case in Figure 7. 

When v = 7/3, the width of the double-spiral is 6ir. This means the first orbit 
starts towards the negative x-axis, winds around one and a half times before reaching 
the positive axis, and then winds one and a half times more to reach the origin again. 
The precession is 57r, so the trajectory closes after two orbits. We show this case in 
Figure 8. However, on this scale it is impossible to see the tight winding of the spirals 
near the origin. Therefore, we demonstrate it with expanded views in Fig. 9. 

As v — > 2 from above, the spirals get tighter and tighter. Their angular widths, 

= 2ir/\v — 2|, become larger and larger as v — > 2. 

6 Unbound Classical Trajectories, v < 2 or \i < 
6.1 The infinite spiral: v = 2 or fi = 

When v reaches 2, there is a singular change. First, the double spiral becomes infinite 
in angular width. But also, the joining of the two sides of the double spiral at p = 1 
and ip = breaks down. It is as if a tightly-wound double spring broke. The ends 
spiral out to infinity, as we discuss below. 

The v = 2 is a special case in other ways. For us, with our E = condition, it 
is the demarcation between the bound orbits, v > 2, and the unbound orbits, v < 2. 
But also, for v = 2 and E ^ 0, it is the demarcation between the E < bound 
solutions and E > unbound solutions. These latter are Cotes' spirals |2||. The 



E 7^ solutions are given in Ref. [25 



Returning to E — 0, Eq. (0) shows that, for fi = or v = 2, there is no natural 
length scale. This is because both the centripetal and the external potentials have 
the same power dependence on r: 



U(r) = (L 2 /2m- 7 )^ . (61) 
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Thus, for all r, the effective potential U(r) is either repulsive (and positive-valued) 
or else it is attractive (and negative- valued). 

For U(r) > there is no solution for E = 0. However, for U < the energy- 



conservation condition (|22|) with E = gives 



= ( 2 m 7 /L 2 - l) r 2 = A V , (62) 

so that 

dr/dip = q=Ar . (63) 
Therefore, the two possible solutions are 

r = r exp [=)=A (y? — ip Q )} , — oo < ip < oo . (64) 

These correspond to orbits which pass through the point (r , ip ) and spiral inwards 
or outwards for the minus or plus signs, respectively. For example, a complete, 
counterclockwise orbit, starting from the initial value r = r in = r(ip in ), going into the 
origin, and then going out of the origin towards the final value r = r/ = r(ifif), is 
given by 



r 

? in 

r 



exp [-A (ip - <f in )) , if in < ip < oo , (65) 



— = exp [+A (ip — iff)] , — oo < ip < iff . (66) 

However, it is important to note that Eqs. (^) and (^) are only special cases of Eq. 
(|64]) . For example, Eq. (|64D, with the plus sign, can also describe an r starting at 
r(tpo) = r which then goes out to infinity as ip — > oo. 

Despite the infinite spirals in the example of Eqs. (|65D and (^), the journey in 
and out takes only a finite time. To see this, consider the time dependence of r. From 
the energy-conservation conditions (|5|) and (|5"2"D, we have 

, 2 _ 2 7 L 2 L 2 [2m^ A \ 1 ( LX\ 2 1 
mr 2 m 2 r 2 m 2 



V L 2 ) r 2 I m ) r 2 ^ ^ 



Therefore, 



and 



—— = 2rr = =F (6^ 

dt m 



r 2 = r 2 T—(t-t ) . (69) 
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This tells us that, in spite of the infinite spiralling, a particle moving inwards reaches 
the origin, r = 0, from any point Tq in a finite time, given by 



At = t ~ to= 2L\ r2 °- (70) 

Finally, substituting Eq. (|64j) into Eq. (|69|) gives the general time dependence of cp 

as 

t ~ to = T^(r 2 - r 2 ) = T ^|{exp[ T 2A(^ - ^)] - 1} . (71) 



Because of our choice L = z, in Eq. ( |7T| ) we are assuming that y> increases monoton- 
ically with time. 

6.2 General unbound orbits: v < 2 or pi < 

When the potential parameter z/ just leaves that of the infinite spiral, that is, when 
one barely has v < 2 or p < 0, there is another change. Although the two ends of the 
entire trajectory still reach to infinity and the spirals in and out almost have infinite 
angular widths, the distance of closest approach jumps from p = to p = 1. 

As the value of v decreases, the value of the angular width of the trajectory, now 
given by = vr/|/x|, also decreases accordingly. By the time v — 1, the angular 
width has decreased to 2tt. Eventually it becomes less than n, meaning the orbit 
comes in and out in the same half plane. That is, when v < 0, the force is repulsive. 
In the next section we give some illustrative examples. 



7 Examples of Classical Unbound Trajectories 

7.1 Kepler-like potentials: 0<^<2or-l<,u<0 

When < v < 2, the repulsive centripetal barrier dominates at small r whereas the 
attractive potential V = —^jr v dominates at large r. A typical shape of the effective 
potential is sketched in Fig. 10. It is familiar from the Kepler problem. Therefore, 
for < v < 2, the E = classical orbits are all unbounded. The distance, a, defined 
in Eq. (|T2| ) now has a completely different interpretation. It is now the distance of 
closest approach. Even so, the formal solution remains valid for negative values 
of p. 

As a first example consider the case v = 3/2 or p = —1/4. This orbit has a 
total angular width of 4tc. It is shown in the two drawings of Fig 11. The large-scale 
first drawing shows the trajectory coming in from the top, performing some gyration, 
and going out at the bottom. The small-scale second drawing shows the trajectory 
winding around twice near the origin, with the distance of closest approach being 
one. 
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A second example is the exact Kepler potential, 1/ = 1 or // = - 1/2. Eq. (f26|) 
gives 

p-V 2 = cosip/2 , (72) 

so that 

- = (cosy>/2) = . 73 

P 2 

This is the famous parabolic orbit for the Kepler problem with E — 0. This orbit is 
shown in the first drawing of Figure 12. The parabola yields an angular width of 2tt, 
as it should. 

7.2 The straight line: is = or [i = — 1 

If we formally set v = in the expression ([|), we get a negative constant potential 
V(r) = —7. Therefore, in this case the force vanishes and we have a free particle. Its 
orbit must be a straight line. However, Eq. (^4|) shows that one still has the same 
type of solution, Eq. fl2~Tp. Here it is 

p = [cos LpY 1 , x = r cos ip = a . (74) 

This is the equation for a vertical straight line that crosses the x-axis at x = a, as 
required by the initial conditions. This orbit is shown in the second drawing Fig. 12, 
it subtending an angular width of 7r from the origin. 



7.3 Repulsive potentials: —00 < v < or —00 < \x < — 1 

For v < the potentials V(r) in Eq. (|j) are repulsive and negative- valued for all 
r > 0, with y(r) going to —00 at large distances. Since both the potential, V(r), 
and the centripetal potential decrease monotonically, the effective potential has no 
minima or maxima. Even so, for E = these unbounded orbits behave qualitatively 
like those for < v < 2. The distance of closest approach again obeys the formula 
(|T2|) and the solutions are given by the same expression (^), which is valid for all 

The solution (|27j ) for the first orbit (<fl = 0) gives 

p = [cos fjup] 1 ^ = [cos . (75) 

This shows that the particle is at r = a when if = and it is at infinity when 
ip = ±7r/(2|/i|). Thus, these orbits become narrower as becomes larger. This is 
similar to the case of the width of the petals of the bound orbits. 

Hyperbolic orbits: The most famous special case of these potentials is the "in- 
verted" harmonic-oscillator potential, with v — \i — — 2. The orbit is given by 
p = [cos2<^] -1 / 2 , so that 

1 = — cos2y? = — (cosV - sinV) = — -. (76) 

cr a z a z a z 
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Thus, the trajectory is a special hyperbolic orbit, whose minor and major axes are 
equal, b 2 = a 2 . (In fact, every v = — 2 solution, i.e. for arbitrary E ^ 0, also yields 
a hyperbolic orbit, but with b ^ a.) We show this orbit as the third drawing in Fig. 
12. Now the angular width has decreased to n/2. 

As the last case, we consider the orbit for v = —4 or /i = —3. This orbit is shown 
in the last drawing of Fig. 12. The orbit subtends an angle of 7r/3, again as it should. 
One sees that as v becomes more and more negative, the orbits will become narrower 
and narrower. This is just as in the bound case, where the petals became narrower 
and narrower as v became more and more positive. 

8 Comments on the Classical Problem 

It is well-known that classical orbits precess for general central potentials. There 
are two exceptions, the Kepler and the isotropic harmonic oscillator potentials |4j]. 
Further, these central potentials are known to have additional conserved quantities: 
the Runge-Lenz vector, A, in the Kepler problem and the quadrupole moment tensor, 
Qij, in the harmonic oscillator problem f26|, p7|. These two examples show that there is 
a close connection between the absence of precession and the existence of dynamically 
conserved quantities. 

On the other hand, it is also well-known that for a specific energy to angular- 
momentum combination, general central potentials can have closed orbits. But these 
combinations have different unique values, depending upon the potential. 

Therefore, it is interesting that the class of E = solutions for the potentials (|T|) 
with v > 2, have a countable number of examples with well-defined, closed trajec- 
tories. They are all the potentials with rational values of v. One of them, the case 
v — 4, has zero precession. 

All the E = with v > 2 solutions are bound and go through the origin. The 
v < 2 solutions are unbound, and go to infinity. (The v — 2 case is special, and does 
both.) 

We have seen that, with the v > 2 singular potentials, the physical solutions differ 
from more familiar solutions in that, at the origin, we must paste together single-orbit 
solutions which have different phase-shifts, ip^, at r = 0. The discontinuity at the 
origin of the phase shifts, given by Eq. fl39|), is fundamental [ES|. It is necessary 



because the bound orbits pass through the origin, where the potential is singular. 
This causes the second order time- derivative of the position, i.e., the acceleration, to 
be infinite at r = when v > 2: 



£ *(*) = * = 1,2,3. (77) 



Therefore, irrespective of what coordinate system we choose, at the origin we must 
still paste together different solutions of Newton's equation which depend on different 
initial conditions, as exhibited in Eq. (j3~ip . 

As we have just observed, when crossing the center of singular potentials, a particle 
will have infinite momentum and kinetic energy. This leads to a violation of the virial 
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theorem, even for these bound trajectories. Formally applying the virial theorem 
29], |30| to power potentials leads to the following well-known equation: 



(T) = ~^(V) , 



(78) 



where (T) and (V) are the time averages of the kinetic and potential energies. The 
above equation leads to two famous relations: (T) = — (V)/2, for Kepler elliptical 
orbits (E < and v = 1), and (T) = (V) for the harmonic oscillator (E > and 



Since (T) > 0, Eq. (|79|) immediately shows that the virial theorem ( f78|) is violated 
by all the -E = solutions of power potentials with v ^ 2. Such a violation of ( f78f) 
is expected for infinite orbits, such as the parabolic orbit of the Kepler problem and 
all other E — 0, v < 2 orbits. What is of more interest is that the bound solutions, 
v > 2, also violate the theorem. In Appendix C we will go into more detail on this 
point 

In summary, the classical E = solutions for the power-law potentials exhibit a 
fascinating set of properties. In paper II we will show that this characteristic is also 
true for the quantum solutions. 





(79) 
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Appendix A: Zero Angular-Momentum Solutions 



Here we discuss the simpler E = L = solutions, where the effective potential 
U (r) is equal to the potential V(r) itself. A particle moves radially, in a straight line, 
between zero and infinity, with no turning points at any finite radius, r ^ 0. 

From the energy-conservation equation (](]) one has that, for E = L = 0, 



dr 
dt 



± 



27 
m 



1/2 



V(t). 



10) 



Unless v = — 2 this is equivalent to 



d(ryl 2+1 _ d (2 1 \ 1 / 2 



(v/2 + l)dt 



rn 



which has the solutions 



r(t) 



v+2 

V ± 



27 
rn 



1/2 



+ l)(t-t c 



2 

I/ + 2 



f ^ -2. 



52) 



where ro = r(to). 

The solution (|82"D does not depend explicitly on the initial velocity v = v(to) 
because r and t>o are not independent. However, because of the condition E = 0, 



using Eq. (|8~CiD yields ± (27/m) 1 / 2 = v r^ 2 . Substituting this relation into Eq. (j32j) 
and expanding yields 



r(t) = r + u (t - t ) - - t ) 2 + 

4 r 



"T2 - #~ to) + 



J3) 



This expansion is absolutely convergent for |t>o(t — to)| < r o- This power series termi- 
nates, if the index v is equal to negative integers, v = 0, —1, —2, • • •. In particular, for 
a constant force, v = — 1, we get the familiar result, r{t) = r + v (t — t Q ) + |a(t — 1 ) 2 , 
where the constant acceleration is given by a = 7/m = Vq/(2to). 

Eq. (^) shows that for v > —2 it takes an infinite time to travel between r 7^ 
and r = 00, but only a finite amount of time to reach the origin r = 0. This finite 
arrival time is given by 



At = t-tr 



r 



\(u/2 + l) v \ 
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Note that At is equal to the time required by a particle moving with an effective 
constant velocity, 

v eff = \{u/2 + l)v \, (85) 

to reach the origin from r = r . In particular, for the free particle, v = 0, the above 
effective speed ( jS5"l) is equal to the actual constant speed v , as it should. 

However, for v < —2, the power of the large square brackets in Eq. (|8~2"D is 
negative. This yields the opposite situation: Now, it takes an infinite time to reach 
r = from any r$ ^ 0. But it takes a finite time, given by fl84|), to travel between 
any tq ^ and infinity. The latter result is consistent with the result of Sec. 4.4, 
that this travel time is finite for L^O. 

Finally, the special case v = — 2 yields the differential equation 



dr 
~dt 



± 



m ) 



(86) 



The solutions are therefore, 



r(t) = r exp 



m J 



r exp 



r / 



57) 



The Taylor expansion agrees with Eq. ([8SD with i/ = —2. In fact, the exponential 
solution (|S7| ) follows from the general solution (|g2"D, by using lim e _ >oc (l + ex) 1//(E = 
exp[x]. Note that for v = —2 the travel time between r 7^ and either r = or 
r = 00 becomes infinite. This is consistent with v = —2 being a transition between 
the regimes where the travel time is infinite to reach the origin {y < —2) or infinity 
(v>-2). 

Therefore, for the L = case, the complete journey from the origin to infinity 
takes an infinite amount of time for all indices —00 < v < 00. 



Appendix B: Similar Orbits and Scaling 

The power potentials have the property that for each solution, ri = ap((p), there 
are infinitely many similar solutions, r\ = \ap((p). The periods, r, the angular 
momenta, L, and the energies, E, of these similar orbits are related by powers of the 
scaling parameter, A, as follows: 



Li 



E 1 



This scaling law can be demonstrated in various ways. For example, Landau and 
Lifshitz use the Lagrangian to show it |14[]. We shall use Newton's equation directly. 
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For the forces corresponding to the power potentials of Eq. (SI) , we have 



d 2 r dV(r) vy „ . - 

m— = -^r = f-r . 89 

dt 2 dr r v+1 v 1 

If we choose an arbitrary scale for the distance, 

p = r/a , (90) 

scale the time as 

T = t/h, *! = a^VWr , (91) 

and multiply Eq. ([89"D by t\/{ma) = a^ +1 /7, we obtain, instead of Eq. (|89|) , a 
dimensionless equation of motion, 

d 2 p f?d 2 f a^ +1 1/7 

-r 



(iT 2 a rft 2 7 r^ +1 



1/ 



^ • ( 92 ) 

Therefore, the solutions of Eq. (R^) must have the following form: 

r = r (t) = a p{ V {T)) = a p{ V {t/t x )) , (93) 

where p(ip) describes the same functional dependence of p on ip for all similar orbits. 

Thus, if we change the length scale by a\ = Xa and the time scale according to 
Eq. © ; 

t x = a 1 : vj2 ^f 1 = X 1+u/2 a 1+u/2 ^Jm/j = X 1+U,2 ti , (94) 

we obtain similar solutions, which differ in size and period. The other scaling ratios 
in Eq. (jS8|) follow simply by dimensional analysis: 

h = ^ 2 h = (A) 2 A-( 1+ ^ = A 1 -/ 2 , etc.. (95) 



The scaling property ( |9"3"| ) is a powerful result. For example, it immediately leads to 
Kepler's third law for the ratio of the periods, t x /ti = A 3 / 2 , and to the independence 
of the harmonic oscillator period on the amplitude of the motion. It also predicts the 
a— dependence of the unit of time, To, that we used in Sec. 4.4. 

Appendix C: Violation of the Virial Theorem 
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Clearly, the violation of the virial theorem by E = 0, v > 2 bound orbits is 
connected to the infinite values of the potential and kinetic energies at r = 0. But to 
better understand the phenomenon, let us return to the derivation of Eq. (|78| ) and 
investigate the relevant assumptions. 

One often starts by formally differentiating the scalar product G = p • r with 
respect to time p9| : 

G = pr + pr = Fr+— pp. (96) 

m 

Applying this result to power potentials, where 

F = -V V(r) = -jur-^r , (97) 

one obtains 

G = vV(r(t)) + 2T(t) . (98) 
Integrating Eq. (|9~8| ) with respect to time, we find 



G(t 2 ) - G{h) = [ 2 G = (t 2 - t x )[u{V) + 2(T)] , (99) 
Jti 

where, in Eq. (|99|), (T) and (V) again denote time averages of T and V, here over 
the time interval (ii,^)- 

Eq. ( p9|) is valid, as long as the particle stays within one orbit (i.e. one petal) 
during the time interval. But Eq. Q9"9| ) is violated when the particle passes through 
the origin. Then G(t) = p • r changes its value from — oo to +oo at the origin, r = 0, 
meaning the derivation of the virial theorem breaks down. 

To demonstrate that p • r becomes infinite at r = 0, we can show that p goes to 
infinity stronger than r goes to zero. From angular-momentum conservation, we have 

G(t) = p ■ r = rpcosa = L cosa/| sina| = L cot a , (100) 

where a is the angle between the vectors p and r. The angle a changes quickly from 
7i to at the r = crossing, since there r reverses its direction from being antiparallel 
to p to being parallel to p. Accordingly, cot a passes from — oo to +oo. 

In a similar way we can extend these results to E < orbits when v > 2. For 
clarity, we summarize our observations: 

The E = orbits of the power potentials V = — r )/r u with v ^ 2 must either go to 
infinity or pass through the origin. In addition, for v > 2, every bounded orbit with 
E < goes through r = 0. 
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Figure Captions 



Figure 1: For v = 4 a plot of the effective potential, U, as a function of r. (See 
Eq. ©.) f/(r) = [l/p 2 -l/p 4 ]. 

Figure 2: A plot of, P^, the precession per orbit, as a function of za 

Figure 3. A plot of the dimensionless orbital period T u = t u /t as a function of v. 

Figure 4. The first three orbits for v — 8. Each orbit is precessed — 27r/3 from the 
previous one, so that by the end of the 3rd orbit, the trajectory closes. In this and 
later figures, we show cartesian coordinates for orientation. 

Figure 5: The first four orbits for v = 6. Each orbit is precessed — 7r/2 from the 
previous one, so that by the end of the 4th orbit, the trajectory closes. 

Figure 6. The orbit for v — 4. It is a circle, and repeats itself continually. 

Figure 7: The first two orbits for v = 3. Each orbit is precessed tt from the pre- 
vious one, so that by the end of the 2nd orbit, the trajectory closes. 

Figure 8: The first two orbits for v = 7/3. Each orbit is precessed 5tt from the 
previous one, so that by the end of the 2nd orbit, the trajectory closes. 

Figure 9: Close-up details of the beginning of the first orbit for v = 7/3. The orbit 
starts at <p = The four graphs show the evolution for (i) — 3tc < (p < — 2.9l7r, 

(ii) — 3tt < if < —5tt/2, (iii) — 37r < (p < —2tt, and iv) — 3n < if < —it. 

Figure 10: The effective potential, U(p) = 1/p 2 — 1/p, as a function of p. 

Figure 11: A large-scale view, and a small-scale view near the origin, of the tra- 
jectory for v = 3/2. 

Figure 12: From left to right, he trajectories for the cases i) u — 1, ii) v = 0, iii) 
v = —2, and iv) v = —4. The curves are labeled by the numbers v. 
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